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Summary. Relativistic outflows represent one of the best-suited tools to probe the 
physics of AGN. Numerical modelling of internal structure of the relativistic outflows 
on parsec scales provides important clues about the conditions and dynamics of the 
material in the immediate vicinity of the central black holes in AGN. We investigate 
possible causes of the structural patterns and regularities observed in the parsec- 
scale jet of the well known quasar 3C 273. We compare the model with the radio 
structure observed in 3C 273 on parsec scales using very long baseline interferometry 
(VLBl) and constrain the basic properties of the flow. Our results show that Kelvin- 
Helmholtz instabilities are the most plausible mechanism to generate the observed 
structures. 



1 Introduction 

3C273 is the brightest quasar known. Its relative proximity {z = 0.158, [1]) 
has made it a paradigmatic object studied throughout the entire spectral 
range [2]. One of the most impressive features of 3C 273 is an apparently one- 
sided relativistic outflow emanating from the quasar nucleus [2] and extending 
up to a deprojected distance of about 60kpc. The jet in 3C273 has been 
observed on parsec scales using ground and space VLBI [3, 4, 5, 6]. Space 
VLBI observations made at 5 GHz with the VSOP^ revealed a double helical 
patterns inside the parsec-scale jet in 3C273 [5]. Linear perturbation analysis 
of Kelvin-Helmholtz (K-H) instability [7] applied to model these structures 
yielded an accurate estimate of the bulk parameters of the flow [5] . A different 
interpretation suggests that a helical magnetic fleld could generate such a 
structure [6]. We study the first possibility. A bulk Lorentz factor W = 2.1 
obtained in [5] is below the values inferred from superluminal motions observed 
in the jet {W = 5 — 10, [4]), which suggests that the K-H instability develops in 
a slower underlying flow, while the fast components are shock waves inside the 
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jet. Four K-H modes were found acting on parsec scales in the jet in 3C 273: the 

second helical body mode with a wavelength of 2 mas, the first elliptical and 
helical body modes, both with wavelengths close to 4 mas, and the elliptical 
surface mode, with a wavelength of 12 mas. Additionally a structure with a 
wavelength of 18 mas was interpreted as an externally driven helical surface 
mode. Our previous works [8, 9] have shown that numerical simulations can 
be used effectively to connect the linear perturbation analysis with studies of 
non-linear regime of K-H instability. Our aim here is to use the jet parameters 
determined in [5] as initial conditions for numerical RHD simulations of a 
perturbed jet and compare the numerical results with the observed structures. 

2 Numerical simulations 
2.1 Simulation 3C 273-A 

We start with a steady jet with a Lorcntz factor W ~ 2.1, a density contrast 
with the external medium rj — 0.023, a sound speed Cgj = 0.53 c in the 
jet and Cs.a = 0.08 c in the external medium, and the perfect gas equation 
of state (with adiabatic exponent 7 = 4/3). Assuming an angle to the line 
of sight 9 = 15°; the observed jet extends up to « 170 pc. With the jet 
radius i?j = 0.8 pc [5], the numerical grid covers 211 i?j (axial) x 8i?j x 8i?j 
(transversal), i.e., 169 pc x 6.4 pc x 6.4 pc. The resolution is 16cells/i?j in the 
transversal direction and 4cells/_Rj in the direction of the flow. A shear layer 
of 2 i?j in width is included in the initial rest mass density and axial velocity 
profiles to keep numerical stability of the initial jet. Elliptical and helical 
modes are induced at the inlet. 

Frequencies of the excited modes are computed from the observed wavelengths 
of modes, A°'^^, corrected for projection effects, relativistic bulk motion and 
wave pattern speed, Vyj = 0.23 c, according to w = 2 it / \*^'^°^ , where 



Three wavelengths have been identified in the simulation: a helical Af ™ = 4 i?j 
perturbation, a helical Al™ = 25 i?j perturbation and an elliptical A3 ™ = 50 i?j 
perturbation. The wave speeds 11^,2 — 0.38 c and 11^,3 ~ 0.2 c are estimated 
for the Al™ and A3™ modes, respectively. An upper limit for v^^i is given 
by the flow speed (0.88 c). With these wave speeds, A|™ would be observed 
as a 2.27 mas structure, Af™ as a 3.37 mas structure and A3™ as a 5.5 mas 
structure. These are similar to the shorter-wavelength structures found in the 
observations (2.27 mas vs 2 mas, 3.37 and 5.5 mas vs 4 mas). It remains to be 
studied why the simulations do not reproduce readily the longer modes. The 
longest mode, with the wave speed given in [5], corresponds to a simulated 
wavelength of 150 -Rj, which is difficult to reproduce even in the grid of this 
simulation, in particular when shorter harmonics grow fast and disrupt the 
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flow. The disruptive evolution may be caused by a number of factors, including 
the absence of a stabilizing magnetic field in the simulations. 

2.2 Simulation 3C 273-B 

In the second simulation, we try to study the efl'ect of precession and injec- 
tion of fast components on the observed structures in the jet. The precession 
frequency is given by the observed 15 yr periodicity of the position angle 
variations in the inner jet ([4]). The frequency of ejections of components is 
set by the reported 1 yr periodicity ([3]). The duration of each ejection is 
estimated to be 2 months, from the approximate inspiralling time from an 
orbit at ~ 6Rg around a 5.5 10® Mq. Velocity of the components is taken as 
constant, with the mean value of those given in [4], i.e., Lorentz factor ~ 5. 
The components are treated as shells of diameter 0.5 i?j ejected at the center 
of the transversal grid. The numerical grid for this simulation covers 30i?j 
(axial) X 6i?j X 6i?j (transversal), i.e., 24 pc x 4.8 pc x 4.8 pc. The resolu- 
tion of the grid is 16 cells/i?j in both transversal directions and 32 cells/iij in 
the direction of the flow. 

We find two structures in the simulation: a pinching perturbation with a 
wavelength of 0.4 i?j caused by the injection of components, and a helical 
perturbation with a wavelength of 4.0 Rj associated with the precession. The 
wave speed associated with these structures is < 0.98 c. This upper limit for 
the wave speed results in observed wavelengths (from Eq. 1) of 0.6 and 6 
mas, much smaller than the 2-4 mas and 18 mas expected from the obser- 
vations. This results further supports the identification of the shorter modes 
with K-H instability modes. We also find that the 15 yr precession period 
cannot account for the observed longer, 18-mas wavelength. Thus, either the 
longer-wavelength structure is driven by a different mechanism, or the 15 yr 
periodicity in the source is not associated with the precession. 
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